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Abstract 

Bulk M0S2, a prototypical layered transition- metal dichalcogenide, is an indirect 
band gap semiconductor. Reducing its size to a monolayer, M0S2 undergoes a transi- 
tion to the direct band semiconductor. We support this experimental observation by 
first principles calculations and show that quantum confinement in layered d-electron 
dichalcogenides results in tuning the electronic structure at the nanoscale. We further 
studied the properties of related TmS2 nanolayers (Tm = W, Nb, Re) and show that 
the isotopological WS2 exhibits similar electronic properties, while NbS2 and ReS2 
remain metallic independent on size. 

1 Introduction 

Layered transition-metal dichalcogenides (LTMDCs) of TmX2 type (Tm = Mo, W, Nb, 
Re, Ti, Ta, etc., X = S, Se, Te, etc.) have been studied extensively on the experimental 
and theoretical basis for the last 40 years.-"— Molybdenum disulfide (M0S2) is a prototypical 
LTMDC, which is composed of two-dimensional S-Mo-S sheets stacked on top of one another, 
as shown in Fig. [TJ Each sheet is trilayered with a Mo atom in the middle that is covalently 
bonded to six S atoms located in the top and bottom layers. The bonding between the 
adjacent S-Mo-S sheets is much weaker and they are held together by van der Waals forces 
leading to the two-dimensional character of M0S2. 

Due to the weak forces between the sheets and the anisotropic character of LTMDCs, 
shearing takes place more easily even under high pressure, leading to lubricant applica- 
tions.-»i^ Other applications, such as catalysis, optoelectronics and photovoltaics, have been 
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Figure 1: The atomic structure of layered transition- metal dichalcogenides of TmX2 type 
(Tm - transition metal, X - chalcogenide) . Different sheets of TmX2 are composed of three 
atomic layers X-Tm-X, where Tm and X are covalently bonded. Sheets are held together by 
weak van der Waals forces. 



proposed and investigated as well.->^i^»^""— A single layer of M0S2 was recently used to real- 
ize a field-effect transistor with Hf02 as a gate insulator.— The transistor exhibits a room- 
temperature current on/off ratio exceeding 1x10* and mobility comparable to the mobility 
of thin silicon films or graphene nanoribbons. Two-dimensional sheets of layered materials 
can be produced nowadays by e.g. liquid exfoliation, what was successfully performed in case 
of transition-metal dichalcogenides by Coleman et al.— 

Electronic structure of various bulk TMDCs and a monolayer of M0S2 has been previ- 
ously studied throughout ab initio calculations using plane waves approach or local basis 
functions.-'^»^>^'^ M0S2 is an indirect band gap material in its bulk form, which recently 
was shown to become a direct band gap semiconductor when thinned to a monolayer.— 
Therefore, few-layer-MoS2 materials can have their band gaps tuned and their values shift 
from the bulk value by various amounts due to quantum size effects. 

In this paper, we have studied size-dependent electronic properties of TmS2 type, where 
Tm = Mo, W, Re and Nb. The first principle calculations were performed using localised 
Gaussian basis functions and compared to the available experimental data and plane wave 
calculations. The results show that the indirect-direct band gap transitions also hold for 
WS2. The band structures of ReS2 and NbS2 reveal metallic character and therefore, no 
transition appears. It can be expected that the size- dependent phenomenon of band gap 
transition will occur also for other dichalcogenides, like MoSe2, WSe2, MoTe2 or WTe2. This 
study is in progress and will be the subject of next communication. 
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Table 1: Calculated and experimental lattice parameters of hexagonal transition-metal 
dichalcogenides in the form of TmX2 (Tm = Mo, W, Nb, Re; X = S, Se). Results ob- 
tained at the DFT/PBE level. In parenthesis, data obtained at the DFT/PBEO level are 
given. 
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M0S2 


3.173 (3.143) 


12.696 (12.583) 


3.160 


12.295 


WS2 


3.164 (3.139) 


12.473 (12.380) 


3.154 


12.362 


NbSa 


3.332 (3.313) 


12.106 (12.074) 


3.310 


11.890 


ReSa 


3.300 (3.275) 


12.724 (12.148) 







2 Methods 



In this work, we have studied LTMDCs of TmS2 type, where Tm = Mo, W, Re and Nb. 
All structures have hexagonal symmetry and belong to the PG^/mmc space group. ReS2 is 
a triclinic system but for the purposes of a direct comparison we have used hexagonal sym- 
metry. The monolayers and polylayers were cut out from the fully optimised bulk structures 
as (0 1) surfaces. Different numbers of layers were considered in this studies: mono-, bi-, 
and quadrilayers, as well as, 6- and 8-layered structures. 

First-principle calculations were performed on the basis of density functional theory 
(DFT) as implemented in the CRYSTAL09 code.— The exchange and correlation terms 
were described using general gradient approximation (GGA) in the scheme of PBE (Perdew- 
Burke-Ernzerhof)— and PBE hybrid (PBEO) functionals.— The following basis sets were 
used: Mo_SC_HAYWSC-311(d31)G_cora_1997 (for Mo atoms), W_cora_1996 (for W atoms), 
Nb_SC_HAYWSC-31(d31)G_dallolio_1996 (for Nb atoms), Re_cora_1991 (for Re atoms), and 
S_86-311G*Jichanot_1993 (for S atoms). 

Optimisation of initial experimental structures was performed using analytical energy 
gradients with respect to atomic coordinates and unit cell parameters within a quasi-Newton 
scheme combined with the BFGS (Broyden-Fletcher-Goldfarb-Shanno) scheme for Hessian 
updating. The optimised lattice parameters for all the studied materials are given in Tabled! 

The shrinking factor for bulk and layered structures was set to 8, what results in the 
corresponding number of 50 and 30 k-points in the irreducible Brillouin zone, respectively. 
The mesh of k-points was obtained according to the scheme proposed by Monkhorst and 
Pack. — Band structures were calculated along the high symmetry points using the following 
path r - M - K - r. 
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3 Results and Discussion 

We have studied electronic properties of TmS2 layered structures with respect to the size 
of the systems. Bulk structures as well as mono- and polylayers were considered. Figs. |2] 
and O show band structures of M0S2 and WS2, respectively, calculated using PBE functional 
and going from bulk to a monolayer. [For band structures results using PBEO functional see 
Supporting Information]. 
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Figure 2: Band structures of bulk M0S2, its monolayer, as well as, polylayers calculated 
at the DFT/PBE level. The horizontal dashed lines indicate the Fermi level. The arrows 
indicate the smallest value of the band gap (direct or indirect) for a given system. The top of 
valence band (blue) and bottom of conduction band (green) are highlighted (Online color). 

The results show that the bulk M0S2 and WS2 are indirect-gap semiconductors. The band 
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Figure 3: Band structures of bulk WS2, its monolayer, as well as, polylayers calculated at the 
DFT/PBE level. The horizontal dashed lines indicate the Fermi level. The arrows indicate 
the smallest value of the band gap (direct or indirect) for a given system. The top of valence 
band (blue) and bottom of conduction band (green) are highlighted (Online color). 

gap originates from transition from the top of valence band situated at F to the bottom of 
conduction band halfway between F and K high symmetry points. The optical direct band 
gap is situated at K point. As the number of layers decreases the indirect band gap increases 
and becomes so high in the monolayer that the material changes into a 2D direct band gap 
semiconductor. At the same time, the direct gap stays almost unchanged (independent of 
size) and close to the value of the direct band gap of a bulk system. These results are 
visualised better in FigJH where the band gap of M0S2 and WS2 are plotted against the 
number of layers. 
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Figure 4: Calculated direct and indirect band gap values of M0S2 and WS2. The horizontal 
solid lines indicate the band gaps of bulk structures. 



There is a significant difference in the band gap values calculated using PBE and PBEO 
functionals. Comparing the results with the available experimental data, one can see that 
the PBE functional performs better. The indirect band gap of bulk M0S2 is 1.23 eV^ and our 
value of 1.2 eV agrees perfectly at the PBE level, while it is overestimated by around 1 eV 
using PBEO hybrid functional. Similarly, the band gap of WS2, with experimental value of 
1.35 eV,^ is 1.3 eV at the PBE level and 2.4 eV using PBEO. The experimental direct band 
gaps of M0S2 and WS2 are 1.74 eV and 1.79 eV, respectively.^ Our DFT/PBE calculations 
give values of 1.8 eV and 1.9 eV for M0S2 and WS2, respectively. Again, the PBEO results 
are by around 1 eV overestimated. 

Li et al.,— using GGA/PBE functional and the plane wave approach, obtained band 
gap of 0.79 eV for the bulk M0S2. For the same system. Matte et al.— have obtained 
1.1 eV and 1.70 eV for the indirect and direct band gap, respectively, using GGA/PBE 
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and DZP basis set. Using Perdew-Wang exchange-correlation functional and self-consistent 
pseudopotential, Arora et al.-- have obtained the indirect band gap of 1.32 eV for the WS2 
bulk. 

Decreasing the number of layers causes progressive shift in the indirect gap up to 1.9 eV 
and 2.1 eV for M0S2 and WS2, respectively. The change in the indirect gap energy is 
significantly larger than that of the direct gap, which increased by around 0.1 eV. Similar 
results were found by Mak at al.^^ for M0S2. 

The band gaps of bulk structures are in the range of near infrared. Reducing the size of 
the material to just a few layers causes the blue shift in absorption features and shifts the 
band gap to the range of the visible light. These materials can be, therefore, interesting in 
optoelectronics. 

These unusual electronic structures of M0S2 and WS2, and the resulting optical properties 
come from the d-electron orbitals that dominate the valence and conduction bands (see Fig. |5] 
for M0S2). The projected density of states (PDOS) of shows that p-states of sulfur atoms 
hybridize with the d-states of the transition metal atoms at the top of valence band and 
the bottom of conduction band. The core states are dominated by the s-orbitals of the 
chalcogenide atom. 

We have also studied the role of quantum confinement in niobium and rhenium disulfides 
and if their band gaps can be tuned by changing the size of the material. The electronic band 
structures of NbS2 and ReS2 (see Fig. [6]) are different from the discussed above, especially 
close to the Fermi level. Although all the studied sulfides are isotopological, the NbS2 and 
ReS2 materials have metallic character independent of the number of layers. 

The metallic character of NbS2 comes from the fact that the 4(i^2 orbital is half-filled and 
results in a band that crosses the Fermi level in several points in the Brillouin zone. The 
band at the Fermi level is derived primarily from these 4rf-orbitals and becomes separated 
from all other states in a monolayer. 

The structure of ReS2 studied here is of the same symmetry as all the other LTMDCs, 
even though the experimentally obtained material is triclinic. ReS2 has one electron per 
metal atom more than M0S2 and has a metallic character similar to NbS2. The band that 
crosses the Fermi level comes from the Srf-orbitals of Re atoms. Unlike NbS2, this band 
in not separated even in the monolayer and the electronic structure in the vicinity of the 
Fermi level is determined by short-range interactions in the sulfur 3p and rhenium 5d band 
complex. 

For both, NbS2 and ReS2, the electronic structure of a bulk can be resemble by a single 
layer of these compounds. 
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Figure 5: Partial density of states of bulk M0S2 and its monolayer calculated as the 
DFT/PBE. The projections of Mo and S atoms are given together with the contributions 
from 4d and 3p orbitals of Mo and S, respectively. The vertical dashed lines indicate the 
Fermi level. (Online color). 

4 Conclusions 

We have studied layered transition-metal dichalcogenides in form of TmS2 by means 
of GGA-DFT (PBE) and a hybrid functional (PBEO). The former was found to match 
experimental results of the bulk phases within less than 0.1 eV, while PBEO overestimates 
the band gap by 1 eV, but gives otherwise very similar band strcutures. Our DFT, at 
the PBE level, support the recent findings of Splendiani et al.— for M0S2: the material 
changes its electronic properties from an indirect semiconductor in the bulk phase to a 
direct semiconductor in the monolayer, an interesting phenomenon with potential application 



8 



NbS^ bulk 



NbS„ quadrilayer 



NbS„ bilayer 



NbS„ monolayer 



0) 2 - 



0) 2 

CD 



- 0.2 




r MKrr MKrr MKrr mkt 

ReS2 bulk ReS2 quadrilayer ReS2 bilayer ReS2 monolayer 




r MKrr MKrr MKrr MKr 



Figure 6: Band structures of bulk NbS2 and ReS2, their monolayer, as well as, polylayers 
calculated at the DFT/PBE level. The horizontal dashed lines indicate the Fermi level. 



for optical devices. We show that WS2 has very similar properties and can be considered 
as alternative material. On the other hand, both ReS2 and NbS2 systems have metallic 
character independent of the number of layers. 

The LTMDCs, if mixed together, have interesting and promising technological potential, 
especially for nanolelectronics and catalysis. Therefore, we are presently working on the 
properties of doped structures of transition-metal dichalcogenides. 
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